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B decays are a subject of active research since they provide useful information on the dynamics 
of strong and electroweak interactions for testing the Standard Model (SM) and models beyond and 
are ideally suited for a critical analysis of CP violation phenomena. Within the standard model, 
there exist certain relations between CP violating rate differences in B decays in the SU{3) limit, 
as for example A{B'^ — > 7r"'"7r~) = — A(_B'' — > n^K~). The goal of this letter is to study the direct 
CP violation asymmetry in a class of processes where there has been recent theoretical progress, as 
for example the B decays into two light pseudoscalars mesons and into a light pseudoscalar and a 
light vector meson. We identify relations between rate asymmetries which are valid in the SU(3) 
limit in the standard model and we compute SU(3) breaking corrections to them, going beyond the 
naive factorization by using the QCD improved factorization model of Beneke et al.. Finally, in 
some processes as for example BR{B~ rj' K~), we claim that one has to add SUSY contributions 
to the Wilson coefflcients. In these cases, we end with a BR depending on three parameters, whose 
values are constrained by the experimental data. 



I. INTRODUCTION 

As it is known, in the Standard Model, the CP vi- 
olation arises solely from the phase in the 3x3 uni- 
tary CKM matrix and any CP violating observable 
is proportional to Im{VijV^V^jVki) , with i ^ k and 
j I. The SU(3) invariant amplitude for B PP 
and B — s- PV decays, in terms of the tree and penguin 
contributions, are, for example. 



A{B° vr+vr-) 



A{B^ TT+K-'' 



= VubV:^T + V,bV:^P, 

= Vu,v:,T+v,,v:,p, 



with T and P the same in the two processes. Even 
there are no simple relations among the branching ra- 
tios of these decays since the CKM factors in the T 
and P amplitudes are different, one has the following 
relations among the CP violating rate differences, [6], 



A?: 



-A 



'■K+K- ^K+TT- 

and similarly for B PV, where 
Afp = r{B PP) - f (B 



-Af 



PP), 



(1) 



(2) 



while the CP asymmetry is defined as 



Acp = 



T{B PP) - T{B PP) 
T{B PP) + f (B ^ PP) 



(3) 



The most important question now is to establish to 
what precision these relations hold within the stan- 
dard model, or equivalently to estimate the correc- 
tions they might receive from different sources, as for 
example the SU(3)breaking effects. 

In this respect, let us discuss the corrections 
in QCD improved factorization method developed 
by Beneke et al. [3, 4] and compare the results 
with the ones obtained within the so called Naive 
Factorization, [1], which will be briefly presented in 
the next section. 

II. NAIVE FACTORIZATION 



The effective weak Hamiltonian for B 

cays,[l]. 



PP de- 



eff 



CiOl + C2OI + J2 C^O^ + Cr^Ov^ + CsgC 



1=3 



h.c. 



(4) 



where Ap = VpbV*^, with p ~ u,c and q = d (for 
AS* = processes) and q = s (for AS = 1 processes), 



in the case when 



0, becomes 
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eff 



Gp 

71 



A„(Ci05' + C20^) 




h.c. 



(5) 



It is expressed in terms of the Wilson coefficients Ci 
(evaluated at the renormalization scale fi = tub), the 
usual tree level left-handed current-current operators, 
the QCD and electroweak penguin operators, the elec- 
tromagnetic and chromomagnetic dipole operators. 



Within naive factorization, one is able to put the 
matrix element of the Hamiltonian in terms of the 
form factors and decay constant of the meson which 
is factorized as 



{P,P2\H,ff\B) = z^Xp (^i-Ci + fpMl - m?)Fo^^^n^2) + (1 2) , (6) 



r 



where Nc = 3 is the number of colors and we introduce with 
the usual combinations of the Wilson coefficients 



ai = Ci + ^Ci+i (for i = odd), < = , < = ^cbK*. (9) 

Ui = Ci + ^Cj-i (for i = even) . (7) and 



3 

In the concrete case of the Bq tt+tt" and Bq ^^(if) = ^^^"^ w 1.2 . (10) 

TT^K~ processes, the matrix elements respectively (wfc — '7i«)('71m + 

are, [1], 

Using the unitarity relations and introducing the no- 
A{Bo -> TT+TT-) = - i^/.Fo^-(m^)(m| - ml) ^^^'^'^'^ 

v'^ _■ v' 

x{via^-vf[a4 + a^o+r„{ae + as)]}, ^^--R^e'^ -| = i?, e '7 , 

A{Bo^Tr+K-) = -i^fKF^^''{ml){ml-ml) a = ai + aw + riae + as) , p = ai+a (11) 

v2 

x{w* ai — vl [ai + aio + (ae + ^s)]}, (8) the amplitudes (8) get the expressions 
I 



A{Bo ^ 7r+^-) = - i^f^Fi^^ml - ml)vt [-Ra e^^^/S + a] , 

A{Bo - TT+K-) = - i^/^Fo^--(m| - ml)v^, [R, e'^^/J + a] , (12) 
which allow us to write down the CP violating amplitude difference 



G 



2 



(l^r-I^P)B^.. = ^f'M^-)\ml-ml)'Aivt)'R,sm^S. 



{\Af-m^^^^ = - ^fl{Frn\ml-ml)H{vl)^R,sln^5^K 



(13) 



I 

where the quantity 5 = Re(/?)Im(a) — Im(/?)Re(a) is Within SU(3) flavor symmetry, when phase space 
the same. differences are neglected, naive factorization yields the 
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relation, [8], 



and Belle, 



J TT 

f2 
JK 



which can be used to test the SM CP violation, or to 

predict one rate difference if the other one is known. 
In these assumptions, the relation between the CP 
asymmetries, 

^CP(.-.-) = -|5|^^CP(.-i^-), (15) 

for the reported CP branching ratios [Babar] 
Br(7r+7r-) = (5.8 ± 0.4 ± 0.3) x 10'^ 

Bt{-k+K-) = (19.4 ±0.6) X 10"^ , 
leads to the following result, 

Acp{tt+tt-) « -2.2 Acp{^+K-) , 

which does not agree with the experimental data that 
are just emerging from Babar, 

Acp{tt+K-) = -0.107± 0.018 



AcriTT+TT-) = 0.55 ± 0.08 ± 0.05. 



III. IMPROVED FACTORIZATION METHOD 



Let us turn to the improved factorization method, 
(IFM), developed by Beneke et al. [3], which gives a 
systematic and modcl-indcpcndcnt calculation of two- 
body hadronic decays, in the heavy-quark limit. The 
factorization formula, presented in the previous sec- 
tion, is applicable, since the nonfactorizable correc- 
tions are included in the parameters, which have 
imaginary parts coming from vertex corrections and 
penguin contributions. In this approach, the Wilson 
coefficients are calculated at the scale /x — mi, us- 
ing next-to-leading order modified scheme, the elec- 
troweak penguin contributions are considered as next- 
to-leading order and there is also a contribution com- 
ing from the hard scattering with the spectator. 

The IFM formula when both Mi and M2 are Ught 

mesons is 



{MiM2\0i\B) = Fo^^^'fM.J dxTij^^icPMAx) + (1 ^ 2) 

+ fBfMjM2 j dz dy dx T" {x, y, z)(j>B{z)(l)M-, {y)4'M2 (x) ■ 
I 



(16) 



where are the leading-twist light-cone distribution 

amplitudes and the integration is over the momentum 
fractions inside the mesons, T/ includes tree diagrams 
plus corrections (hard gluon exchanges and penguins) 
and r/^ expresses the hard gluon exchange with the 
spectator. For T^ = 1 and T/^ = 0, we recover the 



J 



naive factorization. The meson wave functions will be 

an important source of SU(3) breaking. For the light 
mesons, we have twist-2 and twist-3 distribution am- 
plitudes respectively defined in the following bilocal 
operator matrix elements: 



(M(p)|g(^2)7^759(^i)|0) = -/mP^ [\x e'^-P-'^+^P-'^U{x) 

Jo 

{M{p)\q{z2)ij5q{zim = fnm I dx e'^^P-''+^P-'^Up{x) . 

Jo 



(17) 



where /Um is expressed in terms of the quark masses momentum space, the light-cone projector operator 
as hm = 'm\jl{mi + 7712) and x = 1 — x. In the of a light pseudoscalar meson described by both the 
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twist-2 and twist-3 amplitudes is: 

^(M) = ^ {in5ct>{x) - m 75M^)} , (18) 

where p = p ■ ^. 

We notice that in oi, 02, 03, 04, 05, 07, 09, aio, 
where we have {V — A){V±A), only the twist-2 ampli- 
tude is taken, while in ag, ag (the terms proportional 
with r = 2fj,/mb) only the twist-3 amplitude must be 
considered. The operators ae and as are important 
for penguin-dominant B decays. 

The twist-2 distribution amplitude, (/'(x), has the 
following expansion in Gegenbauer polynomials [2] 

<i){x) = 6x(l-x)[l + aiCp/^'(2.T-l) 

+ a2C|/^(2x- 1) + ...], (19) 

with Cl^^iu) = 3u and C^^^{u) = (3/2) (Su^ - 1), and 
is different for tt and K. For tt, the distribution in x 
must be even because the u and quarks have negligi- 
ble masses and their distributions inside the pion are 

I 



ai(MiM2) 
a2(MiM2) 
a3(MiM2) 
<(MiM2) 
a5(MiM2) 
<(MiM2) 
a7(MiM2) 
af(MiM2) 
a9(MiM2) 
afo(MiM2) 



where Cf = {N^ — l)/2Nc, include the vertex, the 
hard gluon exchange with the spectator and the pen- 



symmetric. This dictates = 0. The coefBcient 
is estimated to be 0.25 ±0.15. For K, the u (or d) and 
s quarks inside the kaon are different, leading to an 
asymmetry in the x distribution. So a non-zero value 
for af- is needed and it is estimated to be 0.05, while 
a2 = 0(2- The leading twist distribution amplitudes, 
valid for /i — > 00, are (l>{x) = 6x{l — x) and (pp{x) = 1. 

The B meson distribution between the heavy quark 
and light antiquark is described by 

<t>B=N ^ I (20) 
[a-^^i. — z) + z^y 

where the parameter a is related to the position of 
the maximum of the amplitude and is very small 
a e [0.05 0.1]. However, since the momentum is al- 
most carried by the heavy quark, one may work with 
a strongly peaked function around zq = Xb /ttib ~ 
0.066 ± 0.029, for As = 0.35 ±0.15 GeV. 

In the this approach, [3], the coefficients, 



(21) 



I 

guin contributions, at = nib, defined as 



Nc 



1 + ^{Vm.+H) 



C3 



= Ca 



= C. 



1 + 



1 



47r 
47r 

An 



{Vm, + H) 



l + ^(^M.+^f) 
1 + 



Nr. 



1 - 6 
1 + 



47r 

Cpa. 



47riV, 
(-12 - Vm^ - H) 
C poi. 



pp 



1 - 6 



Cpas 
An 

Cpa. 



" AnNc^'^^' 
{-12 -Vm,-H) 



1 + 



47r 
An 
An 



pP,EW 

9nNc^^'^ ' 



iVM,+H) 



9nNc 



-P. 



.p.EW 



M2,2 ' 
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Vm 



18- 



dxg{x)4>Mix), 



pp 


= Ci 


- + Gm{ 


Sp) 


+ C3 




+ {C4 + Ce)[{nf- 


2)Gm( 


pP,EW 
M,2 


= (Ci+7VeC2) 




pP 

M,3 


= Ci 


- + Gm{ 


Sp) 





-+Gm(0) + Gm(3) 



dx 

X 



(t>M{x), 



3c; 



77 



dx 
Jo ^ 



4>m{x), 



+ Gm(0) + Gm(1) 



+ (Ci + Ce) {nf -2)Gm{0) + Gm{sc) + Gm{1) 



-2Cl 



eff 



r)P,EW 

M,3 



H = 



47r2 fsf. 



-+Gm{sp) 



3c; 



eff 
77 • 



Ml 



Nc msXEFo^'^^'iO) Jo 



dx 



<t>M2 {x) 



dy 
V 



r 



(22) 



where the parameter 2iJiM/mb coincides with r and 
Si = rrii/ml are the mass ratios for the quarks in- 
volved in the penguin diagrams, namely = Sd = 



Ss = Q and Sc = (1.3/4.2)2. 

Putting everything together, the amplitudes (8) get 
the explicit form 



.Gf 



A(B°^7r+7r-) = -i^(m|-m2)Fo^--(m2)/,[K6V„>i(7r7r) 

+ ypby^dKi-^T^) + «io('r7r) + r^(ag(7r7r) + af (tttt)))] , 

+ Vj,bV;Ml(.^K) + aUnK) + rMnK) + a|(7ri^)))], 
where p is summed over u and c, and consequently, the relation (14) turns into, [6], 



n+K 



1 - 0.748< - O.lOQaJ - 0.00137?^ 



1 - 0.748af - O.lOQaf - QMI^H^k 



(23) 



(24) 



(25) 



pointing out the following SU(3) breaking effects: the 
difference in the decay constants and form factors, the 
difference in the a\ and 0:2 coefficients that appear in 
the twist-2 distribution amplitudes (19) and the il^^ 
and Ht^k contributions (defined in (22)). 



relation: 



The decay amplitudes for K^ir^ and i?° 

K~^K~ can be obtained by using the appropriate tran- 

and by changing l/m^As 
Xb, in Hmi M2 • One gets the same expres- 
sion, (25), and thus we have come to the following 



sition form factor Fq' 

to 1/m^ As 



n+K 



(26) 



These example shows that important SU(3) break- 
ing effects arise from the light-cone distributions of 
mesons in addition to those already present in the 
decay constants. These effects can only be estimated 
with large miccrtainty because the parameters af 2 are 
not well determined at present. Using the currently 
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allowed ranges we find, 

Acpin+TT-) « - (S.ltli) Acp{n+K-) , (27) 

which can also be used to test the SM and the IFM 
to some extent. 

However, relations which are independent of a\ 2 
parameters and decay constants, such as (26), are 
more reliable since they do not receive the main SU(3) 
breaking corrections that we have investigated. 

A. B ^ PV Decays 

When the vector meson is factored out, as in 
B° ^ w+p- , B° ^ K+K*- , 



B° ^ K+p- , B° ^ n+K* 



the decay amplitudes can be obtained by replacing the 



Tk factor with = 



nil, J'k' 



0.3 (and similarly 

for Tp), and by removing the penguin terms P^ff^ 
in the expressions for ag and ag (the vector meson is 
described only by a twist-2 distribution amplitude). 
With all these taken into account, we get, [6], 



ruB Ip 
mB. fh 



1 - 1.25(3!^ - O.lSa^ 
1 - 1.25af * - O.lSaf 



(28) 



Using the central values of the ranges 



When the meson that picks up the spectator is a 
vector, as for example in 



K 



the corresponding decay amplitudes can be obtained 
by replacing the form factor F^~*^ with A^^^ and r 
with — r. In this case, the SU(3) breaking is large and 
estimates are unreliable since the analogue of (29), [6], 



rriB fl 



A 



K* + K- 



rriB. 



K 




1 + 110< + 15.5aJ 
1 + UOaf + 15.5af 



(30) 



contains large coefficient of ai in both the numera- 
tor and denominatorm making a prediction for this 
asymmetry impossible within this framework. On the 
other hand, this provides an opportunity to constrain 
(or even to determine) when the ratio in (30) is 
measured. 



K-f]' 



or. 



0, a% 



0.15, 



af* =0.04, af* =0.10 
and taking fp 0.96/x* , we find. 



A^: 

7r+£ 



K+p- 



■0.95 



■0.95 



Ft 



B^TT 



.F 



B,^K 



EL 

FF 



Bs^K 



r}'K decay which has 
CLEO announced its 
-^n'K-) = {&.bt\l± 



Let us turn now to the B~ — > 

become of a real interest after 
large numerical value BR{B~ 
0.9) X 10~^, which could not be explained by the ex- 
istent theoretical models. As improved measurements 
followed, providing even larger values, (80tg" ± 7) x 
lO-'^ (CLEO) and (76.9±3.5±4.4) x 10"^ (BaBar), in- 
clusion of new contributions for accommodating these 
data has quickly become a real theoretical challenge. 



The relevant decay amplitude for B 



T]'K- 



is, 



(29) [4], 
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A{B- rj'K-) 



.Gf , • 
■ % — ^(to 



+ V,,V;, + a\^{X) + rKialiX) + a|(X)))] 



+ 



1 



(a9(y)-ar(y))(l 



,(31) 



where X = rj'K and Y — Krj' , r' = 2rn^,/{mb — 
ms){2ms) and a = f^//f^i- As it can be noticed, 
the coefficients a, are different for the X and Y final 
states, since thciv depend on the twist-2 and twist-3 
wave functions of the M2 meson, except for the hard 
contribution where the wave functions for both Mi 
and M2 arc involved. 

The twist-2 distribution amphtude (j)K{x) has the 
usual expansion in Gegenbauer polynomials, while the 
corresponding twist-3 amplitude, (f>^, is 1. In what it 
concerns the physical states rj and rj', these are mix- 
tures of SU(3)-singlet and octet components 770 and 
?78 and the corresponding decay constants, in the two- 
angle mixing formalism, are given by, [1], f^, = 63.5 
MeV, = 141 MeV, while the relevant form factor 

for the B ^ Tj' transition is F^^^ = 0.137. Even 
the ry' flavor singlet meson has a gluonic content, [10], 
which could bring a contribution to the wave func- 



tion, this is supposed to be small and therefore we 
employ, in the calculation of T/^', P^, 2 and Pfy^^ in 
ai{Y), only the leading twist-2 distribution amplitude 
(pri' = 6xx. Also, since the twist-3 quark-antiquark 
distribution amplitude do not contribute, due to the 
chirality conservation, the penguin parts in ag(l^) and 
a|(F) are missing. 

In IFM, we get for the B~ K~r}' decay the 
numerical value Br{B Kt]') = 3.65 • 10~^ which 
is comparable to other theoretical estimations, [1, 4, 
10], but is only half of the averaged experimental 
data, suggesting that one has to incorporate take into 
account new contributions in order to increase the 
Br{B Krj') numerical values. 

In this respect, we have employed ing the Minimal 
Supersymmctric Standard Model (MSSM) , by adding 
to the effective SM Hamiltonian the following SUSY 
contribution 



expressed in terms of the gluon and photon operators: 



Osg = ^ mbsa^^l + 75)G''''6 , Or^ = rribsa^^l + j5)F^''b . (33) 



The Wilson coefficients are given by, [5] 

GpiVubV:, + KbK* )m? ""^^ ^ 
GpiVubV:, + VcbV*,)ml rub 



2 

g 

7^ {Msusy) - - ^ ... ... ^T/_.T/* W ^LR — M^), (34) 



where 



Go (a;) = — -7 [22 - 20a; - 2a;^ -|- 16a; ln(a;) - a;^ ln(x) + 9 ln(a;)] , 

3(1 — a;)** 

4a; 

Fo(a;) = -i- ri-F4a;-5a;2-h4a;ln(x)-|-2a;=^ln(a;)] . (35) 

9(1 — '- ■' 



I 

In the above relations, x = m|/m|, with rrig be- ing the gluino mass and rriq an average squark mass, 



Insert PSN Here 



8 International Conference on Hadron Physics TROIA '07, Canakkale, Turkey, 30 August - 03 September 2007 



while the factor S^'' = A'"'/m?, where A*"* are the off- 
diagonal terms in the sfermion mass matrices, comes 
from the expansion of the squark propagator in terms 
of S, for A ^ TO?. In principle, the dimensionless 

quantities S^'', measuring the size of flavor changing 
interaction for the sb mixing, are present in all the 
SUSY corrections to the Wilson coefficients and they 
are of four types, depending on the L or R helicity 
of the fermionic partners. However, one finds that 
{cf^^''},=^, for MsuSY = m,- - 500 GeV and 
X ~ 1, do not bring any significant contribution to 
the branching ratio. The situation looks different in 
what it concerns the SUSY Wilson coefficients (34) 
that, when rUg is of order of few hundred GeV, will 

I 



pointing out, besides the (IFM)-value 3.65 x 10~^, the 
SUSY contribution depending on p and tp. 

A detailed analysis of this formula, suggests 
that one should take p £ [0.005,0.01] and ip G 
[— 37r/4, 37r/4], for accommodating the range within 
the two extreme experimental data, B Rexp{BaBar) = 
7x10-^ and BRexp{C LEO) =8x10"^ For p close to 
the lowest limit of its interval, the predicted BR*"*"^'^- 
values lie below the experimental data, while for p 
moving to the central value and — 7r/4 < ip < 7r/4, 
one gets BRtotai e [7 x 10"^ 8 x lO^^]. 
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